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Abstract 

We present a simple model that describes the motion of a single chemotactic cell exposed to a 
traveling wave of the chemoattractant. The model incorporates two types of responses to stimula- 
tion by the chemoattractant, i.e., change in polarity and change in motility of the cell. The periodic 
change in motility is assumed to be induced by the periodic stimulation by the chemoattractant on 
the basis of previous observations. Consequently, net migration of the cell occurs in a particular 
direction with respect to wave propagation, which explains the migration of Dictyostelium cells in 
aggregation processes. The difference between two time delays from the stimulation to the two 
responses and the wave frequency determined by the frequency of the secretion of the chemoat- 
tractant are important parameters that determine the direction of migration and the effective 
interaction between cells in a population. This result explains the dispersed state of a population 
of vegetative cells and cells in preaggregation without the assumption of a chemorepellent, and 
also explains the commencement of the aggregation. The result is extended to a general fact as 
follows: the temporal oscillation of the magnitude of the random motion for gradient-sensing par- 
ticles induces spontaneous movement, even when the particles are exposed to a periodic wave of 
the chemoattractant, which results in the aggregation or dispersion of the particles communicating 
via their attractant. 
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I. INTRODUCTION 



Several types of cells can sense the presence of extracellular signals and migrate in the 
direction of the concentration gradient of these signals. This phenomenon is known as 
chemotaxis, which plays an important role in a variety of biological systems including mam- 
malian neutrophilis fibroblasts {2], microglia js], and cancerous cells ^|. Dictyostelium 
discoideum serves as an ideal model for the study of chemotaxis. Each Dictyostelium cell mi- 
grates toward the gradient of the chemoattractant, cyclic adenosine monophosphate (cAMP) 
5j. It is well known that a population of Dictyostelium cell exhibits collective behavior. 
When bacterial food is available, Dictyostelium cells live as unicellular amoebae. In the ab- 
sence of food, the developmental phase of the life cycle is induced, that is, the cells aggregate 
into a multicellular slug, and form a fruiting body whose spores germinate into amoebae. In 
the aggregation process, the dispersed starved cells periodically emit extracellular pulses of 
the chemoattractant 6-^, a target or rotating wave of the chemoattractant concentration 
is formed, and the cells migrate toward the center of the wave via chemotaxis 

In this article, we present a highly simplified phenomenological model that describes the 
motion of a single chemotactic cell. The analytical results allow us to explain the migration 
of cells toward the wave source in the aggregation process and to give an explanation for the 
aggregation of starved cells and dispersed state of vegetative cells and preaggregating cells 
in a population of Dictyostelium cells. 

Simplified and generalized descriptions of chemotactic migration are desired not only for 
the understanding of natural sciences but also for various industrial applications. Recently, 
numerous phenomena that emerge through natural self-organization processes have inspired 
industrial applications such as nanofabrication [l^ and robotics 



11 



I2I ]. Dictyostelium is a 



typical example of gradient- sensing systems that facilitate self-organization of the popula- 
tion; therefore, the understanding of the mechanism as simply as possible would contribute 
to the development of artificial systems and related technologies. The mechanism of migra- 
tion presented in this study is sufficiently simple to allow application of the mechanism to 
other gradient-sensing systems. One such system is presented in the last section. 
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II. MODEL 



Consider a single cell in one dimension, exposed to a traveling wave of the chemoattractant 
concentration (see Fig. [^a)) represented by a sinusoidal function: 

S{x,t) = Save + SoscSm{kx + Ut) , (l) 

at position x and time t. Here, k and u denote the wave number and the angular frequency, 
respectively. Save and Sosc are constant parameters which satisfy Sosc < Save- Although it 
has been observed that the typical waveform is not a simple sinusoidal form but a sharply 
peaked one for a population of Dictyostelium cells jo], as we discuss later, the use of the 
sinusoid is not essential to the result. The result qualitatively holds as long as the wave is 
monophasic and periodic. 

The model presented in this study is such that the position of a single cell, x{t), is 
governed by the following equation of motion, 

^a;(t) = xit)VS{x{t - Tpoi),t - Tpoi), (2) 

:= Xave+XoscSm[kx{t - Tmot)+(^{t " ^mot)]- (3) 

Here Xave, Xosc, Tmot and Tpoi are constant parameters which satisfy Xosc < Xave- 

This model equation is constructed as follows (Fig. [^b)). We use a widely-accepted as- 
sumption that the velocity of the cell is proportional to the spatial gradient of the chemoat- 



tractant concentration 
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15j. The proportional coefficient, Xi is referred to as chemotactic 
coefficient in what follows. Recent studies on intracellular signal transduction pathways of 
chemotactic cells have shown that the stimulation by extracellular chemoattractant triggers 
three responses which should be separately recognized, namely, directional sensing, polar- 



ization, and change in motility |l6l-ll8| . In terms of those studies, from the viewpoint of the 
modeling of the motion of a cell, we assume that the responses induced by the stimulation 
by a traveling wave consequently can be divided into two responses. One of them is the 
set of the directional sensing and the polarization, that is, the spatial localization of several 
signaling proteins along the plasma membrane because of the presence of the spatial gra- 
dient of the chemoattractant concentration in the environment, which results in pseudopod 
extension to a particular direction. We assume that the factor proportional to VS* describes 
the magnitude of this response. The constant parameter, Xpo/, is introduced as the time 
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FIG. 1: (a) The situation considered in this study. A single chemotactic cell is exposed to a 
chemoattractant wave propagating from right to left. The cell senses the gradient of the chemoat- 
tractant and the temporal change in the concentration, (b) Illustration which schematically shows 
the construction of the model. The stimulation by the traveling wave of the chemoattractant 
induces two responses, change in polarity and motility, with time delays, Tpoi and Tmot- 

required from the stimulation to this response. Another response is the change in motility, 
namely, the magnitude and frequency of the extension of pseudopodia. In the observation 
of cell locomotion, the motility would be recognized as the quantity which indicates how 
fast a cell translocates on the substrate. Thus, we assume that the chemotactic coefficient, 
X, is the factor which quantify the level of this response; therefore, it is assumed that the 
chemotactic coefficient oscillates in response to the oscillation of S. The constant parameter, 
Tmoti is introduced as the time required from the stimulation to this response. 

The assumption of the oscillatory behavior for the chemotactic coefficient is also based 
on remarkable previous studies on cell locomotion. Soil et al. determined the time series of 
the instantaneous velocity of a single Dictyostelium cell in detail when the cell was exposed 
to temporally periodic change in the chemoattractant concentration in the absence of the 
19 -El- that the magnitude of the instantaneous velocity 



spatial gradient 



oscillated at the same frequency as the extracellular chemoattractant concentration. This 
result implies that the temporal change of the chemoattractant concentration induces the 
change in the magnitude of the motility, which means here how fast a cell translocates. When 
a cell is exposed to a periodic traveling wave of the chemoattractant concentration, there is 
also the oscillation of extracellular chemoattractant concentration around a cell. Therefore, 
the model assumes that the chemotactic coefficient oscillates. In addition, because a sig- 
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nificant constant phase difference was observed between the oscillation of the extracellular 
chemoattractant concentration and that of the instantaneous velocity in the experiments, 
the time delay, Tmot, corresponding to the phase difference is introduced in the model. 

The assumption in this model provided above can be described in another way from the 
theoretical point of view. The equation of motion of a single cell is generally given by 

,l. = F + «, (4) 

as long as we neglect the inertial effect of a cell. Here, F denotes external forces, ^ denotes 
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23| . and the coefficient 7 is referred 



forces which induces random locomotion of a cell 
to as the friction coefficient. The equation reads 

= 7-'i^ + 7~'e (5) 

In the absence of the spatial gradient of the chemoattractant and other external forces, the 
governing equation becomes 

I. - 0-'?. (6) 



According to the experiments carried out by Soil et al. 19N21|. the temporally periodic 



and spatially uniform stimulation by the chemoattractant induces the periodic change in 
the instantaneous velocity, ^x, at the same frequency with a constant time delay. Here, we 
attribute the oscillation of the velocity to the oscillation of 7"^. On the other hand, when a 
cell is exposed to a traveling wave of the chemoattractant concentration S, F is proportional 
to VS", that is, 

7|. = VS + «. (7) 

Omitting the random component ^ since we are interested in the time-averaged quantities 
in this study, Eq. @ is reduced to 

|. = XVS. (8) 

where x is proportional to 7^^. Because the magnitude of the extracellular chemoattractant 
concentration periodically changes when the cell is exposed to a periodic traveling wave of 
the chemoattractant, we use the results obtained by Soil et al., that is, we assume that the 
chemotactic coefficient x oscillates at the same frequency as the wave. Incorporating the 



time delays to the responses, we obtain the model equation ([2]) in consequence. Although 
the periodic function that precisely describes the chemotactic coefficient has never been 
determined in the presence of the spatial gradient of the chemoattractant concentration, we 
choose a sinusoidal function for simplicity. As discussed later, this choice of the sinusoidal 
function is not essential to the result; the result holds as long as the chemotactic coefficient 
is periodic function. 

In summary, we postulate in the modeling that a response, change in how large or fre- 
quently pseudopodia are extended, can be distinguished from a response, determination of 
the direction of locomotion; the level of the pseudopod extension is affected by the temporal 
change in the magnitude of the extracellular chemoattractant concentration whereas the 
direction decision is affected by the spatial gradient of the extracellular chemoattractant 
concentration; and the composition of these factors determines the instantaneous velocity 
of a cell exposed to a traveling wave. Explicit incorporation of the two time delays, Tpd 
and Tmot, is an essential feature of this model. It should be noted that this model is not 
another expression of models presented thus far in which the periodic change in motility is 
considered but an essentially different model which provides another result on cell migration. 
The significant difference between them is discussed later. 



Abbreviation 


Meaning 


Used value 


Ref. 


27r/k 


Wave length of the chemoattractant wave 


2000 [fim] 


m 


2it/uj 


Periodic time of the chemoattractant wave 


7 [min] 


m 


kXaveSosc 


Averaged velocity of the single cell 


lOlfim/min] 


[25] 




Amplitude of the oscillation of the velocity of the single cell 


10[fj,m / min] 


[25j 


Tmot 


time delay from the stimulation to the change in motility 


6 [mm] 


[19] 


Tpol 


time delay from the stimulation to the change in polarity 


1 [min] 


* 



TABLE I: Model parameters used in the numerical calculations whose results are shown in Fig. 
2. All of them are typical values for Dictyostelium cells in aggregation process except for Tmot- 
Although Tpoi has not been measured thus far, it is estimated with the model in this study (see the 
Discussion in the main text). 
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III. RESULT 



We immediately see that the propagation of the travehng wave induces net migration of a 
cell to a particular direction in average. We can safely approximate — 7v,oz) = x{t — Tpoi) = 
x{t) because the velocity of the traveling wave, w/A; ~ 300 [/im/mm] |6|, is much greater 
than that of a cell, which varies between and 20 [/im/min] y^]. Thus, the velocity of the 
single cell is governed by 

^a;(t) = kXaveSoscCOS [kx(t) + L0(t - Tpol)] 

+ {kXoscSosc/'^) sin [2kx{t) + u{2t - r^ot - vO] 

-{kXoscSosc/'^) sin [u}{T^ot - vO] • (9) 

On the right hand side of this equation, the third term is a constant whereas the first and 
second terms are periodic with respect to time as long as x{t) is treated as a constant during 
one wave period. Therefore, the third term plays a greater role in the translocation of a cell 
than the first and second terms. The elimination of these periodic terms yields the averaged 
velocity of the cell, v, as 

V = -{kXoscSosc/2) sin [u^Tmot - Tpol)] ■ (lO) 

We see that the cell migrates toward or away from the wave source on average when 
sin [u{Tmot — Tpol)] < or sin [u^Tmot — Tpol)] > 0, respectively. Thus, the frequency of the 
wave and the difference between the two time delays, Tmot — Tpoi, are important parameters 
that determine the direction and the magnitude of the net migration. 

The differential equation ([2]) is numerically solved without any approximation. The result 
is illustrated in Fig. |2l Parameter values used here are summarized in Table [B We see that 
the averaged velocity well describes the net migration of a cell. 

Even if the waveform is not represented by a simple sinusoidal function, u and Tmot — Tpoi 
play important roles in determining the direction of the net migration. Consider a more 
realistic case in which the traveling wave of the chemoattractant concentration is a sharply 
peaked one represented by 

S{x, t) = So exp[Si sin(A;x + out)], (11) 

where 5*0 and Si are constant parameters, and the chemotactic coefficient is given by Eq. Q. 
These functional forms of S{x,t) and x{t) describe the result of the experiment conducted 
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by Soil et al. more precisely 19l-l2l|. In the same manner as the case of the sinusoidal 
traveling wave, it is proved by means of the Fourier expansion that 



V = -{XoscSqSiH /-k) Sm.[uj{Tmot - Tpol)], (12) 

where / := j_^dt^T^ e-xY>{Sit) > 0. We see that the direction of net the migration of 
a cell is determined by the sign of sm[ijj{Tmot — Tpoi)]] therefore, the the frequency of the 
traveling wave and the difference of the two time delays are important parameters again. 

More generally, let us discuss the reason for the emergence of the net migration with 
referring to Fig. [3l If the shape of the wave of the chemoattractant concentration S is 
monophasic and periodic, its spatial gradient VS becomes a biphasic wave whose time 
average is 0. If the time delays are such that the peak response of the oscillating motility 




FIG. 2: The time course of the velocity v{t) and the displacement of a cell x{t) obtained by 
numerically solving the differential equation ([2]) without any approximation. The thin curved line 
at the bottom S{t) represents the chemoattractant concentration, which is shown as a reference in 
order to be compared with the phase in the oscillation of the velocity and the displacement. The 
straight line represents the approximate displacement vt calculated with the averaged velocity v in 

Eq. (nni). 
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x{t) appears when V5'(t — Tpoi) is positive, the time average of the product of x(^) and 
VS{t — Tpoi), which is proportional to v, becomes positive. This imphes that the cells 
migrate toward the source of the traveling wave (Fig. Conversely, if the time delays are 
such that the peak response of the motility x(^) appears when VS(t — Tpoi) is negative, the 
time average of the product of x(t) VS{t — Tpoi) becomes negative. This implies that 
cells migrate in the direction of the wave propagation (Fig. [3]d). Thus, even if the wave is 
not represented by a simple sinusoidal function, the difference between the two time delays 
from the stimulation to the two responses determines the direction of the cell migration of 
the cells. Moreover, as long as we fix the two time delays, Tmot and Tpoi, Fig. and Fig. 

illustrate situations in which the frequency of the traveling wave, u, is small and large, 
respectively. This indicates that the direction to which a cell migrates also depends on the 
wave frequency u. Thus, the mechanism of directional migration of a cell is fundamentally 
the same as long as a wave is monophasic and periodic. The frequency of the traveling wave 
and the difference of the two time delays are important parameters which determine the 
direction of the net migration. 

Note that, if there are no such time delays, that is, if Tmot = Tpoi = 0, the time- average of 
xVS* becomes independently of the functional form of S. The directional net migration 
of cells has been successfully derived using several models in which time dependence of 



15|, 



26|-|29|. On the basis of the discussion 



chemotactic coefficient is taken into account 
provided above, we can conclude that all these models have been constructed such that the 
time delays are implicitly incorporated in order to introduce an appropriate phase difference 
between the oscillation of the stimulation and the responses. 



IV. DISCUSSION 

This result ( ITOi) and ( |T2l) are consistent with the aggregation of Dictyostelium cells in 
that the each averaged velocity v is positive for typical parameter values in the aggregation 
process shown in Table [H which indicates that cells migrate toward the center of a target or 
spiral wave of cAMP concentration. 

The numerical result shown in Fig. [2] well duplicates previous observations of the instan- 
taneous velocity of a single Dictyostelium cell during the aggregation processes 19l-|21| 



m 



that a period in which a cell migrates rapidly alternates with a period in which a cell is 
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FIG. 3: (a) The relationship among the extracellular chemoattractant concentration S, its gradient 
with a time delay VS'(t — Tpd), the chemotactic coefficient with a time delay > 0, and the 
velocity v{t) oc xi't)'^ S {t — Tpoi) , in the case which results in net migration toward the source of the 
traveling wave, (b) The case which results in net migration away from the source of the traveling 
wave. The Figure (a) and the Figure (b) also illustrate the situations in which the frequency of the 
traveling wave is small and large, respectively. Distinction of these situations is important not only 
to see the direction of the net migration of cells but also to infer the effective interaction between 
cells (see the Discussion for detail). 

virtually stationary. In particular, the result successfully reproduces the observation more 
precisely than another model |l5| in that the time courses obtained in the experiments 
displays complex behavior which cannot be constructed by a simple sinusoid such as the 
characteristic small peaks in the stationary period. 

The result provided above implies that there exists an effective interaction between cells, 
because, in the case of Dictyostelium, each cell secretes the chemoattractant, the secreted 



chemoattractant propagates in space, and it serves as a stimulation to nearby cells jsj. More 
precisely, if the frequency of the secretion of the chemoattractant, which determines the 
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wave frequency, and the time delays satisfy the condition that a cell migrates toward or away 
from the wave source, there is effective attraction or repulsion between cells, respectively. 
This fact indicates that, for systems in which cells communicate with each other using 
the chemoattractant as an autoinducer, the frequency of the secretion and the time delays 
determine the effective interaction between the cells. In aggregation process, because the 
parameter values satisfy v > 0, the effective interaction between cells is attraction. This 
attraction is interpreted as the microscopic reason for the aggregation of cells. 

Increase of the frequency of the secretion of the chemoattractant during aggregation has 
been observed s], The averaged velocity flTUl) and ffT^ indicate that, as the frequency 
of the secretion of the chemoattractant, u, increases, the sign of v changes from negative to 
positive at a critical frequency Uc, which implies that the effective interaction between cells 
changes from repulsion to attraction. Even if the wave form is not so simple as Eq. ([1]) or Eq. 
(fTTj) . we can infer that the transition occurs according to the discussion provided above and 
Fig. [HI Therefore, it is possible that the dispersed state of cells before the commencement 
of aggregation is due to the effective repulsion between cells induced by the lower frequency 
of the secretion of the chemoattractant. 

In terms of the transition between the dispersion and aggregation in a population of cells, 
we can estimate Tpoi for Dictyostelium cells. From Eq. (fTOj) and Eq. (fT2l) . the frequency 
of the secretion of the chemoattractant at which the interaction between cells changes from 
repulsion to attraction, cUc, satisfies Tmot — Tpoi = tt/wc. According to a previous experiment, 
the time period of the periodic secretion at the commencement of the aggregation was 
27t/uc ~ 10 [min] [s^. Then, since Tmot ~ 6 [min] 19], we obtain Tpoi ~ 1 [min]. Although 
this parameter Tpoi has not been directly measured in the presence of temporally periodic 
spatial gradient of the chemoattractant concentration thus far, the result is consistent with 
an estimation based on the time which requires for the localization of proteins on the plasma 
membrane 3l|]- Inversely, this consistency suggests the existence of the transition at around 
2-71 /uc ~ 10 [min]. 

Repulsion between Dictyostelium cells in a nutrient environment was reported several 



decades ago 



331] . Although experimental studies on vegetative cells are less than starved 



cells in the developmental phase, it is certain that vegetative cells are single-celled amoeba 
and never commence the aggregation. The present model provides us an explanation for 
this state. It has been observed that right after the nutrient deprivation (within 5 hours 
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precisely), cells pulse chemoattractant randomly every 15 to 30 min 8|. It can be inferred 
that vegetative cells secrete the cAMP at the same or lower frequency. According to the 
discussion above and in terms of Fig. |31 the interaction between cells is repulsion at those 
frequencies. Note that, even if the wave consists of a single pulse, there arise also repulsion. 
Thus, it is possible that the effective repulsion mediated by the chemoattractant allows 
vegetative cells to live as unicellular amoebae and never start to aggregate. 

The fact that the interaction between cells may become repulsion indicates the possibility 
of the non-existence of the chemorep client. Repulsion between vegetative Dictyostelium cells 



was originally named negative chemotaxis [33']. Negative chemotaxis was later defined as a 
variant of chemotaxis in which cells migrate away from higher concentrations of a chemical 
source. Such a chemical substance is called chemorepellent. Due to this historical back- 
ground, it was believed that vegetative cells secrete a chemorepellent and repel each other 
owing to the chemorepellent. A potential candidate for the chemorepellent was reported 
recently {mI- However, on the basis of the present model, such a chemorepellent is not 
required to explain the repulsion between cells. Even if there is no chemorepellent, and 
cells exhibit only positive chemotaxis, that is, x is always positive, attraction and repulsion 
can be switched depending on the frequency of the secretion of the chemoattractant. We 
should note that this discussion just gives an possible scenario which explains for behavior 
of a population of vegetative cells because it has been commonly recognized and accepted 
that the intracellular molecular dynamics such as signal transduction and gene expression is 
largely different between vegetative cells and starved cells. The response of vegetative cells 
to cAMP should be investigated more in the future in order to confirm the discussion here. 

It seems that the model presented here is just an another expression of what has been 
presented thus far in terms of time delay or phase difference of oscillations. However, there 
is an essential difference between them. Although it has been also assumed that the motility 
of cells under a traveling wave of the chemoattractant changes, the motility is believed to in- 
crease during periods in which the concentration of the external chemoattractant increases, 
and decrease during periods in which the concentration of the external chemoattractant 



decreases 



19|. This has been considered as the reason why the traveling- wave chemotaxis 



occurs. However, on the basis of this conventional idea, a cell exposed to a traveling wave 
always migrate toward the wave source independently of the frequency of the wave; therefore 
the effective interaction between cells becomes always attraction according to the discussion 
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presented above. This does not explain for the transition from the dispersed state to the 
aggregation of cells. It should be experimentally confirmed in the future that to which direc- 
tion cells migrate when cells are exposed to a traveling wave whose frequency is much lower 
than the wave emerging in aggregation processes. Note that in such an experiment the wave 
is not necessarily a periodic wave if the durations between peaks are long sufficiently. If those 
cells migrate away from the wave source, the model presented in this study would be valid. 
Or this model will be supported if it is observed that cells right before the commencement 
of aggregation significantly migrate in the direction of the wave propagation. 

V. CONCLUDING REMARKS 

A phcnomcnological model which describes the motion of a single chemotactic cell is 
presented and analyzed in this study. The result of the analysis implies that systems in which 
cells communicate with each other using the chemoattractant as an autoinducer exhibit 
aggregation and dispersion of cells. The frequency of the secretion of the chemoattractant 
and the time delays from the stimulation by the chemoattractant to the responses of polarity 
and motility determine the interaction between cells, attraction or repulsion. 

It must be controversial whether or not the responses to the stimulation by the chemoat- 
tractant can be incorporated as simply as this model because the intracellular signal trans- 
duction consists of actually complicated processes. However, the incorporation allows us 1) 
to understand vegetative cells living as unicellular amoebae and not starting aggregation 
in terms of repulsion between cells, 2) to understand commencement of the aggregation of 
starved cells in terms of the transition from the repulsion to attraction, 3) to derive the cell 
migration toward the source of the target or spiral wave of the chemoattractant in aggre- 
gation, 4) to describe the time course of the velocity of a cell in aggregation more precisely 
that other models presented thus far, and 5) to estimate the time which requires for the 
polarization. 

The model in this study just provides at present a possible mechanism which is used in 
cell migration and induces dispersion and aggregation of cells. The vafidity of this model 
should be confirmed by some experiments in the future. 1) r^ot and r^o/ are regarded as 
constant in this study. It should be confirmed whether or not values of r„iot and Tpoi depend 
on the frequency of the stimulation, and taking the dependence into consideration also leads 
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to the same results. Note that, even if those parameters depend on the frequency, the 
mechanism proposed in this study works as long as there remains the critical frequency at 
which the transition occurs. 2) We have used the value of Tmot that was measured in the 
absence of spatial gradient of cAMP in this study. Whether or not the change in T^ot is 
the same in the presence of the spatial gradient of cAMP. 3) Because it has been reported 
that many of responses occur within a few minutes 16| , why the change in motility requires 
such a long time, about 6 [mm], should be elucidated in terms of the intracellular dynamics. 
4) Net migration of cells exposed to a traveling wave whose frequency is small should be 
investigated in order to compare with models presented thus far (see the preceding section 
for detail). As long as we refer to the results obtained in this study, the intriguing question 
arising next is how cells regulate the frequency of the secretion of chemoattractant depending 
on the environment, that is, how they choose sufficiently low frequency in order to live as 
unicellular amoebae when they can survive and grow thanks to enough food, and how they 
choose sufficiently high frequency in order to start development when they cannot survive 
alone because of starvation. The answer to this question will elucidate the mechanism of 
the adaptation of cells to their environment and the strategy taken for the species survival. 

From the viewpoint of physics, it should be remarked that the random motion of a cell 
plays an important role in the generation of directional migration even though the ran- 
dom motion itself is nondirectional, and does not seem to contribute to tactic behavior. In 
addition, it should be emphasized that the composition of symmetric conditions induces 
asymmetric behavior, that is, the composition of the periodic stimulation and the nondi- 
rectional fluctuation of the position generates directional net migration. This spontaneous 
symmetry breaking leads to aggregation or dispersion in a population. 

Let us further generalize the results. Consider an isolated element that moves in response 
to the spatial gradient of a field. When such an element is exposed to the temporally periodic 
change of spatial gradient of the field and the motility oscillates in the same frequency, spon- 
taneous directional motion occurs on average. The direction of the motion is determined by 
the frequency of the change in the gradient and the phase difference between the oscillation 
of the motility and that of the response to the spatial gradient (Fig. Hj). Moreover, a group 
of elements, each of which emits a field with the same frequency (e.g., quantum dots exposed 
to electromagnetic waves JssI]), is expected to exhibit aggregation or dispersion as seen in 
Dictyostelium cells. In the future, such a simple mechanism can potentially facilitate the 
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V{t) ~ sm{ujt) 




71 









T{t) - sin(tj(t - r)) 



FIG. 4: An example in which directional motion occurs with the same mechanism described by the 
proposed model. Consider an electrically charged particle in a solution whose temperature can be 
controlled. The particle is sufficiently small to move randomly in terms of the thermal fluctuation. 
If the particle is exposed to an AC electrical field and if the temperature is changed with the same 
period as the AC field with a constant time delay r, the particle moves in a particular direction 
on average. The time delay r determines the direction of movement. 

development of self-organizing systems in which elements communicate or interact with each 
other via a field. 
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